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Pauli exclusion principle can lead to a supression of transport through double quantum dots, even if both 
quantum dots are connected by a resonant driving ac field. In this work, we study how this effect, known as 
spin blockade, can be removed by means of the influence of the driving field on the contacts to the leads. 
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1 Introduction 

One of the main properties of few electron quantum dots is the supression of charge current when the 
energy cost of introducing an extra electron in the system is large due to the charge repulsion inside the 
quantum dot (QD). In such a case, the system presents Coulomb blockade, showing characteristic curves 
where current only flows for certain values of the chemical potentials] 1 ]. Between these peaks, current is 
allowed only through second order processes, which we do not treat here. 

However, by introducing an external AC field, the electrons obtain (from the interaction with the field) 
enough energy to fulfill the energy requirements and tunnel through the system [ 2] . In this way, a finite 
current may appear even in the zero bias configuration (pumping regime). These photon-assisted tunneling 
(PAT) processes through the contact barriers have been studied in AC driven single quantum dots 1 3 1 but are 
usually neglected in the theoretical study of double quantum dot (DQD) resonant pumps|4], where only 
interdot PAT is considered. 

If one takes into account the spin of the electron, another interesting effect takes part in two-site systems 
like, for example, a DQD having one level each. If an electron is trapped in one of the quantum dots, 
transport is only available through the two-electron singlet state of that QD, when an electron with the 
opposite spin tunnels from the other QD and afterwards to the electrode. But, if the other QD is occupied 
by a trapped electron with the same spin, Pauli exclusion principle does not allow the formation of the 
doubly occupied state and, therefore, the current is blocked. This phenomenon is known as spin blockade 
(SB)0. 

In this work, we study a concrete case in which spin blockade may appear in AC driven DQD spin 
pumps[6] and show that PAT processes through the contacts can be important, allowing the trapped spins 
to tunnel out of the system breaking the SB effect. 
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2 Theoretical model 

Our system, consisting in two quantum dots weakly connected in series to two unbiased electron reservoirs 
by tunnel barriers, is described by the Hamiltonian H = Hl + Hr + Hl^r + -H/eads + Ht- We consider 
one level in each QD containing up to two electrons: Hj={L,R} = J2a £ j<y^ja^'j<y + ^j^iT^jJ.' wnere jtj 
is the creation operator of an electron with spin a in dot j and energy Ej a and Uj is the charging energy of 
each dot. The energies Ej a include the Zeeman splitting, Aj, due to the introduction of a magnetic field in 
order to break the spin degeneracy of the different levels, in such a way that we can consider the spin-up 
state as the ground state. Thus: Sj± — £jf + Aj. The reservoirs are described by the term: Hi ea ds = 
S/ e {L R}ka £ ikdl ka dika, where the operator d lka creates an electron with moment k and spin a in lead I. 
Each QD is coupled to the other and to the leads through the terms: Hl^r = —^lr J2a ^La^Ra + h.c. 
and Ht = R}ka(^ c ^ikcr^ la ' + h.c). The constant that describes the tunneling through the contacts, 

7, is asumed to be small and, for simplicity, similar for both QDs and will be treated as a perturbative 
parameter. 

We also introduce an external AC field acting on the energy levels of the quantum dots such that (con- 
sidering h = e = 1): e l(R)o — > £ l(R)o (*) = s l{R)cj ± ^f-costut, where Vac an d w are the amplitude and 
frequency of the field, respectively. The time dependent field will serve as a driver of electrons from one 
QD to the other, allowing the formation of doubly occupied states, which contribute to the current through 
the device. 



2. 1 Master equation 

We study the electron dynamics of the DQD system using the reduced density matrix operator, p = tr^x, 
obtained by tracing all the reservoir states in the density operator of the whole system, % The Liou- 
ville equation, p(t) = —i[H(t), p(t)] gives us the time evolution of the system. Assuming Markov and 
Born approximations [7 1, we derive the master equation for the density matrix elements||8), Pm'mif) = 
(m'\p(t)\m): 

P(t)m'm — ~i^m'mPm'm(t) ~ i[HL^ft(t), p(t)] m i m (1) 
^ ^7n / nPnn{t) ^ F?unPmm(t) I ^m'-m ^m' 'mPvn' 'ra(^)(l ^m'm)i 

in the particle basis: |1) = |0,0), |2) = | t,0), |3> = | J.,0), |4) = |0,T), |5> - |0, |), |6) = | t,T>. 
1 7) = I U>. |8> = I U), |9) - I |,T), |10) = I Tl,0), |H) = |0,U>. |12> = I U,T), |13> = I TU). 
1 14) = I T>Ti}> 1 15) = I I, Tl). 1 16) = I TJ-jTI)' Here, w m > m is the energy difference between the 
states |m') and \m) of the isolated DQD and ft m i m describes the decoherence of the DQD states due to the 
interaction with the reservoir. The time dependence of the energy levels has been transferred to the interdot 
coupling: 



00 

m\H' L ^ R (t)\n) = Ma)e iuut (m\H L ^ R (t)\n), (2) 



v— — oo 



where J„(a) is the ^-th order Bessel function of the first kind, being a = Vac the dimensionless AC 
field intensity. The tunneling rates through the contacts, T mn , are affected by the AC field: 

oc 

T m n = J„ ( -X J Cmf)(WmTi + w), (3) 



v— — oc 



where 



= {f(s)S Nfn>Nn+1 + (1 - /(-£))<W,JV„-i}r (4) 
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are the usual tunneling rates obtained when PAT through the contacts is neglected. Nk — (M^j<r\k) = 
N 3 k is the number of electrons in the system in state \k), /(e) = 1/(1 + e( £ ~ M ^) is the Fermi distri- 
bution function, where j3 — l/kgT and /i is the chemical potential of the leads, and T = 2ir\j\ 2 . These 
transition rates (Eq.©) are related to the decoherence through the relation: JJf2 m / m = §QIfc^ m ' ^km' + 

We obtain the current that flows through the right contact with the relation Ir = Yl m m ' {^m'mPmm — 
r mm ' p m 'm')$N R -i n r , — En^ 1 The s P m currents, Irj and Zr,j,, only account for the processes 
involving the tunneling of spins with up and down polarization, respectively. 



3 Numerical results 

If both reservoirs have the same chemical potential, /i, i.e., there is no bias voltage applied to the DQD, 
and p < Ui + £/ + A/, the non-driven system will be in a stable state if it contains one electron in each 
dot. Additionally, if p > Ur + er, the spin down electron will be the only one able to tunnel out to 
the right lead from the doubly occupied singlet state and the spin up will be trapped in the right QD. 
Thus, the breaking of the spin degeneracy by the introduction of a magnetic field leads to an asymmetry 
in the transport properties of the system, so the spin components of the current, Jr. j and Ir.i, will behave 
differently. 




Fig. 1 Schematic diagram of the non-driven device showing the chemical potentials associated to the transitions 
involving the extraction through the contact barriers of one electron from the doubly occupied state in each QD. Since 
H > Ur +er, the transitions extracting an electron with spin-up polarization from states with two electrons in the right 
dot through the right contact are energetically unavailable, unless they are mediated by the absorption of photons. In 
our configuration, the Zeeman splitting and the energy of the levels are the same in both dots, that is: Al = Ar = A z 
and el = £r- 

Introducing an AC field in resonance with the states | |, f) and |0, Tl) (and, since Al = Ar, also 
with |0, tl) and | 1, J.}), the spin down electron will be delocalized between both quatum dots and there 
will be a finite probability for it to leave the DQD to the right lead. If PAT processes through the contacts 
are not considered, one should expect a net spin down current through the system (through the pumping 
cycle: J., f) |0, Tl) -» {|0, 1)or\ |, tl}} -> I I, T)) but, since the empty left QD can be also filled 
with a spin up electron, the system asymptotically evolves to the state | |, |) (i.e., pa.e(t —> oo) = 1 and 
pij(t — > oo) = 0, otherwise) that leads to SB|6|. 

However, the rates ® allow the "trapped" spins in the DQD to absorb a certain number of photons 
and tunnel out to the leads giving a finite ocuppation probability (through the sequences: | f, T) — > |0, | 
) — * I I) T) or T> T) — ► I T j 0) — ^ | ti I)) to the states | |, f ) and | f, |) which are in resonance with the 
state |0, f |), that contributes to the pumping of a spin down electron to the right lead (Fig. |2). Then, PAT 
through the contacts creates a finite current through the system (Fig. 0, removing the SB. 



© 2003 WILEY- VCH Verlag GmbH & Co. KGaA. Weinheim 



6R. Sanchez, G. Platero, R. Aguado, and E. Cota: Removing spin blockade by photon-assisted tunneling in double quantum dots. 



1 




11 

14 



6 

8 
9 



4 



1 0,4 - 





0, 



5 10 15 20 



time 



Fig. 2 Time evolution (normalized to t — 2n/Q,, being fi = 2Ji(a)tLR the Rabi frequency of the derealization 
processes) of the diagonal density matrix elements ^4,4, p@ t 6, ps,s, P9,9, Pii,n, Pi4,i4 and 1015,15, which describe the 
occupation probability of the states that contribute to the current. We consider the initial condition: pa,a(t = 0) = 1. 
In the case where PAT through the contacts is not considered, we would obtain P6,e(t) = 1, for all times|6|. We do not 
show the occupation probabilities of other states that are not directly involved in the pumping processes. Parameters 
(in meV): t L R = 0.005, T = 0.001, Ul = 1.6, Ur = 1.3, Al = A R = 0.2 (corresponding to a magnetic field 
B « 8T), e L = e R = 0.5, p = 1.9, lo = uu, 8 = Vac/2. 

4 Conclusions 

We have shown that the interaction with the driving field applied on a DQD system affects not only the 
interdot tunneling but also the tunneling through the contact barriers. In concrete, PAT through the contact 
barriers affects the ocupation of the states, giving a finite probability to states that are energetically unstable 
but open new channels to the electronic transport through the system. Therefore, they should be taken into 
account when studying properties of driven quantum dot devices. 
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Fig. 3 Pumped current (normalized to the tunneling probability, T) as a function of the frequency of the AC field (in 
meV). Considering the blocking state, | |, |), if the electron in the left QD absorbs a photon and tunnels to the left 
lead, net spin down current to the right lead is created due to the sequence: | f, f) — > |0, T) — > | J., T) & |0, T J) —> 
{|0, 1)or\ T, Tl)} — * I Tj T)> while spin up current flows in the opposite direction (from right to left) by the sequence: 
| t, T) ->■ |0, t> ->■ II. T) |0, TJ.) & | T, 4) -> I T. TJ.) ->• I T, T) (in this c y de . s P in down current through the right 
contact is also produced). The double arrow (<4>) represents the resonant derealization processes inside the DQD at 
u) ~ Ur = 1.3. On the other hand, if the spin up in the right QD is extracted from | f, f), there is a positive contribution 
to spin up current through the sequence: | |, T> -> I T, 0) -» | T, J) <^ |0, TJ) -» {| T, TJ>or|0, T)} -»■ | T, T>- The 
contribution of this sequence is smaller (since it compites with the sequence | f, |) — > | f , 0) —* \ j, J.) — * | |, J.) — > 
Tj TJ) ~^ I T) T) mat recovers rapidly the state | T, T) without contributing to the current) and is only apreciable at 
high enough field intensities. However, it may be the responsible of the supression of negative spin up current near 
resonance, giving a small bump. Note also that, since so many states are contributing to the dynamics of the system, 
the behaviour of the resonance peaks differs from the typical Lorentzian shape. The parameters are the same as in Fig. 
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